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ALEXEI TSYGVINTSEV 

Abstract. In this work we study the quahtative properties of real analytic bounded 
maps defined in the infinite complex strip. The main tool is approximation by continued 
g-fractions of Wall [12] . As an application, the ABC-flow system is considered which is 
essential to the origin of the solar magnetic field [T] . 



1. Introduction 



O 
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^ ' In 1948 Hubert Wall introduced the particular class of functional continued fractions 

(D '. 

["T I . called gf-fractions. The objective of the present study is to broaden our understanding of 
■ Wall's ideas in the dynamical system theory. 

In this section we will remind the reader of some key facts from the analytic theory of 
Q ■ continued fractions. Let HI = C_ U C+ U (—1, +oo) where C+ = {z E C : Im{z) > 0}, 

C- = {zeC : Im{z) < 0}. 

For arbitrary real sequence Qi E [0,1], i > 1 we call g- fraction the continued fraction 



(N ■ g{z) = {91,92,... \z} = - — , (1.1) 

> : 1+1+1 + 1 

a\ ■ ^ 

\ converging uniformly on compact sets of EI to an analytic function (see [12] ). The map 

■ 9{z) is rational if and only if 9^ G {0, 1}, for some k>l. 
In particular, if = p G {0,1), i> 1 then 9{z) is algebraic and is given explicitly by 

{p,p,---|4 = -— . ^ , ^GC_UC+U(-l/4p(l-p),+oo). (1.2) 

1 - 2p+ a/1 + Ap{l-p)z 

^1 It is known that some ratios of hypergeometric functions can be expressed with help of 

■ (7-fractions. Let a, h, c are real constants satisfying —l<a<c,0<h<Cy^O and 
F{a, h, c, z) be the hypergeometric function of Gauss. Then, as shown in [B]: 

F{a + l,h,c,-z) 
F[a, b, c, —z) 

where 

a+k h+k-1 

Q2k = ; , Q2k-i = ; , k > 1 . 1-4 

^ c + 2A;-l' ^ c + 2A;-2' " ^ ^ 
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In the simplest case: 
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F(l,l,2,-2;) 



ln(l + z) 



(1.5) 



We define the truncated continued (/-fraction as the ra-order approximation of (11. ip : 

(1.6) 



r I T 1 fl-l^; (1 - 9l)92Z (1 - 9n-l)9nZ 

{91,92, ■■■,9n\z} = - — 

1 + 1 + 1 ... 1 



which is a rational function of z analytic in H. 
The next result is due to Gragg [5]: 

Theorem 1.1. Let z G (— l,+oo), then the n-order truncation error satisfies: 

1 



\{9i,92, •••k} - {9i,92, ■■■,9n\z}\ < 
and does not depend on values of 9i, i > 1. 



1 



l + z 



1 - 



1 + 



The real a priori bounds for the (/-fraction are given by the next result: 

Theorem 1.2. ([9J) 

a) Let k = 2n + 1, n = 0,1, then 

Akiz) < 9{z) < Bkiz), -l<z<+oo, 

where 

Ak{z) = {9i,92,...,9k\z}, Bk{z) = {9i,92,...,9k,l\z} . 

b) Let k = 2n, n = 1,2, then 

At{z) < 9{z) < B+iz), 0<z<+oo, 
A-{z)<9{z)<B,{z), -l<z<0, 

where 

Atiz) = {9i,92,-,9k,l\z}, Bl = {9i,92,...,9k\z} , 
and A,=Bt, B' = At . 



(1.7) 



(1.8) 
(1.9) 

(1.10) 
(1.11) 

(1.12) 



Using the above formulas we write below the rational a priori bounds for the (/-fraction 
(11. ip corresponding to k = 1, 2, 3: 
Case k = 1. 



Aiiz) 



1 

1 + 9iz' 



Bi{z) 



l + {l-9i)z 
l + z 



(1.13) 



Case k = 2. 
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At{z) 



Case k = ?). 



'1 - gig2)z + 1 



A2 = , Bc) 



92O- - 9i)z + 1 
{.91 - 9i92 + 92)z + 1 



;i.i4) 
;i.i5) 



Mz) 



B.iz) 



{93 + 92 -9392 - 929i)z + 1 



5'i5'3(l - 92)z'^ + (5-3 + 92 + 91 - 9392 - 9i92)z + 1 

92{l - 5'3)(1 - 9i)z'^ + {l + g2 - 9392 - 9i92)z + 1 
(1 + z){{gi +92 - 9392 - 9ig2)z + 1) 



1.16) 
;i.l7) 




X 



Figure 1. Bounds r{x)Ai{x), r{x)Bi{x) (bold line) and r(x)A^(x) 
r{x)B2{x) (dashed line) for x G (—0.9,6), gi = 0.7, g2 = 0.3, r{x) = 



VT+x. 

The interesting link between (^-fractions and probability theory was reported by Gerl 
One considers the nearest-neighbour random walks X„, n = 0,1, 2,... on Nq = {0,1, 2,...} 
with the one-step transition probabilities pi^k = Prob[X„+i = k \ X„ = i] defined by 



Po,i = 1, Pjj^i = 9j, Pj,j+i = ^-9j 



(1.18) 



with < < 1, j > 1 and Pj^k = in any other case. 

We introduce = Prob[X2„ = | Xq = 0]- the probability of return to in 2n steps. 
The generating function for this sequence 



G,{z) = Y^{^IYp'^,z\ 



:i.i9) 



n=0 
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can be written then as a ^f-fraction: 

Goiz) = {gi,g2, •••k} • 



;i.20) 



In the next section we will examine the relaltion between ^f-fractions and real analytic 
bounded functions. 

2. Functions bounded in the complex strip 

We denote Am,b the set of functions f{z) satisfying the following conditions 

a) f{z) is holomorphic in the infinite strip Sb = {z & C : \Im{z)\ < B} , B > 0. 

b) /(M) C R. 

c) 1/(2)1 < M, M >0,y z e Sb. 



4 fmfs} 
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Figure 2. Map defined by f{z) G Am,b- 
Our present aim is to obtain ^f-fraction representation for the class Am,b ■ 
Theorem 2.1. Let f{z) G Am^b- Then for some /iq > and & [0,1] , k > 1 one has 



f{z) = M 1 



/ioexp(fl) {^i,c/2,.-|exp(f ) -1} + 1_ 
Proof. We introduce the complex domains 

Bm = {z e C : \z\ < M}, M+ = {2 G C : Re{z) > 0} , 
and the conformal maps defined by: 



M-z' 



and 



l(z) = -log(l + 2), / :H^S£ 
V{z) = r\z) = exp (^^) - 1. 



(2.1) 

(2.2) 

(2.3) 

(2.4) 
(2.5) 
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We note that rj is a bijection between M and (—1, +00). 

One verifies that the composition F = m o f o I is holomorphic in H and F{M) C 1H[_|_ 
with F{z) G M for z > —1. Thus, according to theorem of Wall [12], p. 279 F can be 
written as follows 

Fiz) = ^^oVTT-z f ^ , (2.6) 
Jq 1 + zu 

for some nondecreasing real bounded function fj,{u), u G (0, 1) and /iq > 0. 
For / = mT^ o F o 7] one obtains the following formula 



The integral in fl2.7p can be transformed to the continued (^-fraction form 



that together with (12. 7p implies f l2.ip . Proof is finished. □ 
Let 

= j^fiBzM . (2.9) 
To calculate the coefficients Qp in (12. ip one has formulas: 

gp = CMO),0'iO),...,e^^\0)), p>l, (2.10) 

with rational functions Cp determined by calculation of derivatives of both sides of (12.11) 
at z = 0. The recurrent formulas for all Cp can be derived from p^, p. 203. 
Introducing 

0„ = 0W(o) = __/W(o), n>0, (2.11) 
we provide below explicit formulas for /xq, gi, 92- 

fo = ii|^, (2.12) 

_ 1 (16^^ - 8^2 - ^0 + el + el - 89200 - 1)(1 - Op) 

2 ^i-ei + Ae^)iAe,-i + ei) ' ^ ■ ' 

Definition 2.1. We denote by A^^''^ C Am^b, k = 1,2, ... the set of functions for which 
the g -fraction representation (12. ip satisfies the condition 

g.^{0,l}, yt = l,...,k. (2.15) 



6 Alexei Tsygvintsev 

Our aim is to estimate the time of return of f{z) G Am.b to the initial value /(O) i.e 
to study the points Zq G M* such that f{zo) = /(O). To do this we will use the a priori 
bounds fll.l3p applied to the (/-fraction in formula fl2.ip . For p = 2k + 1 one obtains: 



1 



fioVlTv(^Ap{r]{z)) + 1 



<M<,i 

- M - 



, (2.16) 



for z eR. 

\{ p = 2k then 



/io^rT^A+(r/(z)) + l^ 
for z G (0, +oo), and 



<M<n 

- M - 



(2.17) 



/iov/l + ^(^)^p(^(^)) + l, 



for z G (— oo, 0]. 



<M<n_ 

- M - 



In the next theorem, for a given / G ^mbi describe a neighborhood of origin 

in which 2: = is the only solution of /(z) = /(O). 

Theorem 2.2. Let f{z) G a£^^^, /'(O) ^ w/iere 5(1 zs defined by f l^TT]) and f l2T^ as 

function ofM,B, f{0), /'(O). Lei r G M* 6e t/ie j^ozni swc/i that /(r) = /(O). T/ien 



fl0^1 + 7]{z)Bp{7]{z)) + l 



(2.18) 





log 


(^) 


r > — 


TT 







> 



(2.19) 



Proof. One considers fl2.16p with p = 1. We have Ai{0) = -Bi(O) = 1 and define ti, t2 as 
non-zero solutions of the following algebraic equations 



VT+hAi{ti) = 1, VT+t2Bi{t2) = 1, ti, t2 G (-1, +00) . (2.20) 
Simple algebraic calculations show that the only solutions satisfying f l2.20p are given 

1-2(71 ^ 2gi-l 

ti — ■ 

which are related by 



by 



9i 



2 ' ^2 



:i-^?i 



l2 ' 



(2.21) 



1 1 

h^T2~' 



(2.22) 
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Since r]{z) is a bijection between M and (— l,+oo) there exist unique real numbers Ti, 
T2 G M satisfying the following equations: 

viT,)=tu viT2)=t2. (2.23) 

As easily seen from fl2.5p : Ti = — T2 and the proof of fl2.19p follows straightforwardly from 
the formula ([23D. □ 

The next result shows that f{z) G A^-*^, under some conditions on derivatives /*-^^(0), 
p = 0, 1, 2, always returns to the initial value /(O) i.e admits the oscillatory property. 

('2) ________ 

Theorem 2.3. Let f{z) G Ay^, /'(O) 7^ where gi,g2 are defined by formulas f l2.11l) 
and f l27[3|) . fimj) . 

/. VFe assume that the point {gi, g2) G (0, 1)^ belongs to one of the four regions E, F , G, 
H defined by: 

E = {{gi, g2) e (0, 1)^ : > 0, < (71 < 1/2, Q < g2 < 1/2} , (2.24) 

F = {{gu 92) e (0, 1)^ : > 0, < < 1/2, 1/2 < (72 < 1} , (2.25) 

G = {((?!, g2) G (0, 1)' : /^i > 0, 1/2 < (71 < 1, < (72 < 1/2} , (2.26) 

H = {((71, g2) G (0, 1)2 : /^i > 0, 1/2 < (71 < 1, 1/2 < (72 < 1} , (2.27) 

= (1 - ^i)2 - 4(72^2(1 - g2), D2 = gl - 4(1 - g,)\l - (72)^2 . (2.28) 

Let 



c = 


-log{t?) 

71 


> 


if 


{91,92) e E 


(2.29) 


c = 


-log(*f) 

71 


< 


if 


(91,92) e F 


(2.30) 


c = 


-log(t«) 

TT 


< 


if 


(5-1,5-2) e G 


(2.31) 


c = 


-log(tS^^) 

TT 


> 


if 


(51,52) e iJ 


(2.32) 



where t^j \ i,j = 1,2 are defined as functions 0/51,52 by 



,(1) _ 1 -5'! ± V^i ,{2) _ 5i±V^2 



251(1-52) ' '^'^-2(1-51)52- ^'-^^^ 
r/ien there exists r G M* sitc/i i/iai 

fir) = /(O) . (2.34) 

and 

r G (0, C) if C > and r G (C, 0) if C < . (2.35) 
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//. Let e eR* be such that f{e) = f{0), then 



and 



^e(_oo, — log(t«)]U[— log(t«),+oo), zf /'(0)>0, 
n 71 



^G(-oo, — log(tf^)]U[— log(tf),+oo), /'(0)<0. 



TT 



TT 



(2.36) 
(2.37) 



32 







— 








D2 = ^ Y 


/ 
/ 




I — = 


^/ 





0.2 «.4 Oil O.S 

91 



Figure 3. Four domains E, F, G, H in the parameter space {gi,g2) G [0, 1]^. 



Proof. We consider f l2.17p with p = 2 and define the following real algebraic equations 

a/1 + xAj{x) = 1 , xG(0, +CX)), (Ai), 
VTT^B+{x) = l, a;G(0,+oo), (Bi) , 



Vl + xB^{x) = 1, xG(-1,0), (Ai) 
Vl + a;A^(x) = 1, xG(-1,0). (Bi) 



where = 5^, ^2 = Aj. 
Making the change of variables 



x = -l + t, te 



(2.38) 



after some elementary transformations, it is easy to show that equations (^i), (-Bi) are 
equivalent respectively to quadratic equations (^2) and {B2) given below 

Pi{t) = gi{l - g2)t^ - {1 - gi)t + 9192 = 0, teR, (A2) 

P2it)=92{l-9i)t^~9it + {l-9i){l-92) = 0, teR. (B2) 
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Remark 2.1. We notice that P2{t) is obtained from -Pi(t) by transformation 

g,^l-g,, 1 = 1,2. (2.39) 

The polynomial Pi (t) = has two real roots t^i^ , t^2'' ^ 

(1) _ l-g,-y/Dl (1) _ l-gi + y/D; (i) < (i) .2 40) 

25'i(l-5'2) 2gi{l-g2) 
if and only if the following condition holds 

D, = {1- g,f - Ag'Ml - ^72) > . (2.41) 

P2{t) = has two real solutions t^i \t^'' G M 

(2) _ 91 - (2) _ 91 + /D^ (2) , (2) .^x 

2(l-5'i)5'2 2(1-5(1)5(2 

if and only if 

D2 = 9I - 4(1 - 9i)\l - 92)92 > . (2.43) 

Applying the Vieta's formulas to polynomials A2 and B2, and taking into account that 
gi G (0, 1), z = 1, 2 one checks that: 

tf>0, I, J = 1,2. (2.44) 

Case A. Let /'(O) > 0(<(=^ gi < 1/2). Then /(z) is increasing function in the inter- 
val (— e, e) for some small e > 0. We assume that inequality -D2 > holds, so both 
roots t\ and t2 are real. One has -P2(l) = 1 ~ > 0, so, in view of fl2.44p . either 

< tf^ < t^P < 1 (a) or 1 < tf^ < tP (b). One verifies with help of fl2:42|) that (a) is 
equivalent to L2 = gi — 2(1 — 51)52 < and (b) to L2 > 0. Thus, in view of fl2.38p . if 

(2) ^ 

(b) holds, the equation (i?i) will have solution x = — 1 + 1]^ £ (0, +00) and if (a) holds, 
(Bi) will have solution a; = -1 + 1^^^^ G (-1, 0) in view of fl238|) . 

One verifies directly that the condition < 51 < 1/2 implies Di > 0. So, the the both 
roots t^i^ and t^2^ are real distinct numbers. Since -Pi(l) = 251 — 1 < we have < t^i^ < 

1 < t^2^. So, the equation (Ai) will have the unique real solution 1/2 = — 1+^2^^ G (0, +00) 
and {Ai) will have the unique real solution yi = —1 + t\ ' G (— 1, 0). 

Case B. Let /'(O) < 0(-v^ 51 > 1/2). Then f{z) is decreasing function in the inter- 
val (— e, e) for some small e > 0. We assume that inequality -Di > holds, so both 
roots t^i^ and t2^^ are real. One has -Pi(l) = 251 — 1 > 0, so, in view of fl2.44p . either 
< tj^^ < t^^^ < 1 (c) or 1 < tS^^ < t^2^ (d). One verifies with help of (1^^ that (c) is 
equivalent to Li = 1—51 — 251(1—52) < and (d) to L2 = 1—51 — 251(1—52) > 0. Thus, in 
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view of f l2.38p . if (d) holds, the equation (Ai) will have solution x = —l + t\ G (0, +oo) 
and if (c) holds, (Ai) will have solution x = —1 + ^ G (—1, 0) in view of f l2.38p . 



One verifies directly that the condition gi > 1/2 implies D2 > 0. So, the the both 
roots t\ and ^ are real distinct positive numbers. Since -P2(l) = 1 ^ '^Qi < we 
have < tf^ < 1 < t^2^. So, the equation (Bi) will have the unique real solution given 
by ?/2 = — 1 + 4 ^ (0, +00) and will have the unique real solution defined by 
l/i = -l + tf)'G(-l,0). 

Since ri{z) is a bijection of M and (— l,+oo), there exists unique real number G M 
satisfying equation r]{Q = x with x G (— l,+oo) defined above. Then, as follows from 
(EUD, (EUD, there exists r satisfying flCTj) if one of the cases ^Ml-fgMl holds. One 
has r G (0, C) if C > and r G (C, 0) if C < 0. 

Using ?/i 2 defined above, we define 2:1 < and 2:2 > as unique real solutions of 
ri{zi) = yi,i = 1, 2. Let now 6* G M* be such that f{9) = /(O), then 6 G (-00, zi]U[2;2, +00) 
that shows (I2.36p . fl2.37p and finishes the proof. 

□ 

3. Applications to solutions of the ABC-flow equations 

The ABC-fiow is a system of three ordinary differential equations 

/I • < n dx2 dx3 , . 

— — = /Ismx3 + C C0SX2, — — = -D smxi + /lcosx3, — — = C sma;2 + -D cosxi , (3.1) 
dz dz dz 

depending on three arbitrary real positive constants {A, B, C) 7^ (0, 0, 0). This vector field 
appears as an exact solution of the Euler equation without forcing (see [1] for details). It 
is essential to the origin of magnetic fields in large astrophysical bodies like the Earth, 
the Sun and galaxies. 
We define 

6 = max{A + C,B + A,C + B}>0. (3.2) 

Since the vector field of (13. ip is a bounded one it is complete in and hence all its 
real solutions {xi{z),X2{z),X3{z)) are defined for ^ G M. Moreover, the following result 
holds 

Lemma 3.1. Let e > be an arbitrary positive real number. Then every real solution 
{xi{z),X2{z),X3{z)) of (13. ip is analytic function of z in the complex infinite strip 

S, = {zeC : \Im{z)\<b}, 6=1-^. (3.3) 
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Proof. For a given e > and a real triplet [Xi, X2, X^) G M'^, we introduce the complex 
domain T^Xi,X2,X3 

= {(xi,X2,X3) G : \xi - Xi\ < e, i = 1,2,3}. Writing the system 
(13.1 p in the form = fi{xi, X2, X3), z = 1, 2, 3 one easy verifies that 

|/,(P)| <5ch(e), VP=(xi,X2,X3)GP3,^,^^,^3, l<z<3, (3.4) 

in view of the following elementary property: 

Proposition 3.1. The trigonometric functions cosx, sin a; are analytic and bounded in 
absolute value by ch(e) in the complex disk |a; — a;o| < e with a center xq G M. 

We can apply now the classical theorem of Picard (see for example [3]) from the analytic 
theory of ordinary differential equations: 

Theorem 3.1. Let Qj{qi,q2, . . . ,qn) , 1 < j < n be analytic functions in the complex 
domain 

h -Qil <qi, l<i<n, (3.5) 

for some g^' > 0, G C, 1 < z < n. We assume there exist positive constants Q'j > 0, 
1 < J < such that 

\Qj{qi,q2,...,qn)\ <Q'j, l<j<n, (3.6) 

if the conditions (13. 5p hold. 

Then the system of n ordinary equations 

^ = Qj(9i,?2, • • • ,gn), l<j<n, (3.7) 

admits the unique solution analytic in the complex disc 

O-.UKT, r...{|4,...4}, (3.3) 

satisfying the initial conditions qi{0) = qi, 1 < i < n. 
Moreover, 

\qi{z) - qi\ < q'i, zeD, l<i<n. (3.9) 

According to this theorem the solution of (13. ip {xi{z),X2{z),X3{z)) defined by the 

initial condition {xi{zq) , X2{zo) , x^^zq)) = {Xi, X2, X3), zq will be analytic in the disk 

Dzo T' '■ \z — zo\ < T' oi the complex time plane with T' = . Since the system (13.11) 

ch(e) 

is autonomous and 5*;, = ^zoeRDzQ^T' that finishes the proof. □ 

The solutions of the ABC-fiow are defined on the torus = {(xi, X2, x^) : x-i mod 27r}. 
Thus, instead of Xi{z) it is natural to consider the projections sm{xi{z)) : [0, 2ti] — )• [0, 1]. 
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Lemma 3.2. Let {xi{z),X2{z),X3{z)) be any real solution of the ABC-flow fl3.ip . Then 
all functions (f)i{z) = sin{xi{z)), i = 1,2,3 are analytic in Sb defined by (13. Sp and bounded 
in absolute value by M = ch(e). 

Proof. Let Xj{z) G {xi{z),X2{z),X3{z)). Then, (pj{z) = sm{xj{z)) is analytic in Sb as a 
composition of analytic maps. Let Zq G M, then for the complex disc D^gM '■ \z — Zq\ < b 
we have obviously -Dzo.fe C Sb- According to (13.91) W z G Dzo,b '■ \xj{z) — Xj{zo)\ < e. That 
completes the proof in view of Sb = UzoeRD^^^^b and Proposition 13.11 

□ 

The function e(e) = e/ch(e) reaches the unique maximal value for e G M+ which we 
denote Cmax- The direct computation gives Cmax = ^{^max) = 0.6627 for emax = 1.1997. 

We consider an arbitrary solution F^^^^ of (13. ip starting from the plane Xi = 0, and 
defined by initial conditions of the form 

xi(0)=0, X2{0) = a, x3(0) = /3, a, /3 G [0, 27r] . (3.10) 

Let f{z) = (f)i{z) = sin(a;i(2;)). We have /(O) = 0, / G Ajv/,b and the formulas (12. lip . 
dZISD, dlHD give 



1 e 1 
^0 = 0, 01 = — -^(Asin/3 + Ccosa), 6^ = ^ -j—AB cos /3 , (3.11) 
vro ch (e) vr^o^ ch '"^ 



and 



Since / G AM,b, according to Theorem 12. II we have gi^2 G [0, 1]. Indeed, one easy verifies 

(2) 

that the strict inequalities always hold: < gi^2 < 1) so / G Ay^. 

Let r 7^ be such that /(r) = sin(xi(r)) = 0. According to Theorem 12.21 we have the 
following lower bound 



ii_„,Unh(^,, 

7r5ch(e) \7i5ch\e) 



|r| > _ ^ ^ arctanh ( — ^ ^ ^ {Asm(3 + C cos a) ] , (3.13) 



which holds for any (a, /3) G [0, 27r]^ and arbitrary e > 0, 5 is defined by (13. 2p . This means 
that the solution Ta^p starting from the plane Xi = at 2; = can not cross any of the 
planes xi = 0, ivr earlier than permitted by (I3.13p . In practice, one can use a freedom in 
the chose of e in order to make (13.131) optimal. The more precise lower bounds, involving 
g2 i.e the second derivative of / at 0, are given by (12.360 . (I2.37p . 
Now we will analyze the upper bounds given by Theorem 12.31 
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The conditions fl2.24p . f l2.27p . defining the domaines E and F in the parameter space 
(a, /3) G [0,27r]^, are respectively equivalent to 



(the similar descriptions can be derived for the domaines F and G). One can show that 
for any positive constants A,B,C and e the above sets of parameter values [a, (3) are non 
empty in [0, 27r]^ . Indeed, it is sufficient to put « = 7r/2 in fl3.1ip and consider the values 
(/3 — )■ vr, /3 < vr) and (/3 2tt, f3 < 2n) respectively to satisfy (E) and (H). 

The inequalities fl3.14p define those values of for which we can guarantee the exis- 
tence of r > such that sin(xi(r)) = 0. From the dynamical point of view it corresponds 
to returning of the trajectory r^^/j either to the initial plane xi = or to its intersection 
with planes Xi = ±7r. 

The corresponding upper bounds for r can be calculated with help of formulas fl2.24l) . 



We consider the particular case of the ABC-How (13. ip defined by the parameter values 
{A,B,C) = (1 ■ 10-^ 12 ■ 10-^3 ■ 10-^). We take the particular solution r^,/? given by 
(I3.10p with {a, (3) = {tt/2,it — 0.0007). The corresponding parameter values are: e = 1, 
5 = 0.15 ■ 10~^, b = 4320.3618, (91,92) = (0.499987, .494117). 

The lower bound (I3.13P writes as follows: r > 0.137267 and the interval in (I2.36P is 
given by 6* G (-00, -74.8181] U [10.0929, +00). One checks that (5-1, 5-2) e E together with 
flOHj) given by r G (0, 15.2680). 

Performing the numerical integration of the differential system (13. ip one can calculate 
the smallest r* > such that sin(a;i(r*)) = sin(xi(0)) = which corresponds to the value 
r* = 11.6662. Thus, the upper and lower bounds for r* given by Theorem 12.31 are quite 
satisfactory. 

The aim of our study was to elaborate on the ideas of Poincare and Sundman [2] 
from the Celestial Mechanics providing converging time series solutions for the n-body 
problem. Unfortunately, these solutions, though converging for all values of time, have 
very slow convergence. One would have to sum up milliards of terms to gain any significant 
qualitative information above the motion of particles. The present study was designed 
to test the hypothesis that this gap can be overcome by replacing the power series with 
functional continued fractions. We notice that this important issue has not yet been 



{Oi > 0, 02 < 




(3.14) 



(12:2711 and fl233|) . 
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addressed fully in the literature and a number of aspects of the g^-fractions approach 
presented here require further investigation. 
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